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1. Introduction 

Ultracold atoms in 3D optical lattices provide an intriguing environment to study 
strongly correlated condensed matter systems and quantum information. Unique fea- 
tures of these atomic many body systems include the complete control of system param- 
eters, and - in particular contrast to solid state physics - weak couplings to dissipativc 
environments. This so-called quantum lattice gas 012] is described by a Bose or Fermi 
Hubbard Hamiltonian. The high control over the atoms opens the possibility to engineer 
a wide class of interesting many body quantum states. Seminal experiments have already 
demonstrated the superfluid-Mott insulator transition [3], the realisation of ID quantum 
liquids with atomic gases (see also [13 El)? and a Bose spin glass [S]. Here we review 
another recent experiment [5] where we have observed a novel kind of bound state of two 
atoms which is based on repulsive interactions between the particles. These repulsively 
bound pairs exhibit long lifetimes, even under conditions when they collide with one an- 
other. Stable repulsively bound objects should be viewed as a general phenomenon and 
their existence will be ubiquitous in cold atoms lattice physics. Although the experiment 
described here is based on bosonic 87 Rb atoms, other composites with fermions |10| or 
Bose-Fermi mixtures should exist in an analogous manner. Furthermore, repulsively 
bound objects could also be formed with more than two particles. 
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In the following we will first explain the theoretical background of repulsively atom 
pairs. Afterwards we will present the experiments which demonstrate several key proper- 
ties of the pairs. Finally we give a short discussion of how these repulsively bound pairs 
relate to bound states in some other physical systems. 

2. Repulsively bound pairs 

Optical lattices are generated by pairs of counterpropagating laser beams, where the 
resulting standing wave intensity pattern forms a periodic array of microtraps for the cold 
atoms, with period given by half the wavelength of the light, A/2. This periodicity of the 
potential gives rise to a bandstructure for the atom dynamics with Bloch bands separated 
by band gaps, which can be controlled via the laser parameters and configuration, as 
shown in Fig. ^ The dynamics of an atomic Bose-Einstein condensate loaded into the 
lowest band of a sufficiently deep optical lattice E] is well described by a single band 
Bose Hubbard model ^2] with Hamiltonian 

(1) H = -J^blk + | £ % X ) + E 

(ij) j i 

Here bi (b\) are destruction (creation) operators for the bosonic atoms at site i. J/h and 
U denote respectively the tunnelling rate of atoms between neighbouring sites, and the 
collisional energy shift from interactions between atoms on the same site. The resulting 
width of the Bloch band is 4 J, and this single band model is valid because the kinetic 
energy and interaction energy in this system are much smaller than the separation of the 
Bloch bands u>. 

The Bose- Hubbard Hamiltonian predicts the existence of stable repulsively bound 
atom pairs. These are most intuitively understood in the limit of strong repulsive interac- 
tion U >> J (where U > but this energy is still smaller than the separation to the first 
excited Bloch band, U <C ui). If a state is prepared with two atoms occupying a single 
site, \2i) = (S| 2 |vac))/\/2 , then it will have a potential energy offset « U with respect 
to states where the atoms are separated [see Fig. EJleft)]. This state will be unable to 
decay by converting the potential energy into kinetic energy, as the Bloch band provides 
a maximum kinetic energy for two atoms both at the edge of the Brillouin zone given by 
8 J << U. Instead, the atoms will remain together, and tunnel through the lattice as a 
bound composite object - a repulsively bound pair. 

We can observe this nature of repulsive binding in the experiment [see Fig. fright)]. 
After production of the atom pairs, as is discussed in detail in section [S] we allow the 
atoms to tunnel through the lattice along one dimension. If the on-site interaction of the 
atoms is tuned to zero with the help of a Fcshbach resonance, the pairs break up within 
a few ms, corresponding to the tunneling timescale. However, if the effective interaction 
between the atoms is repulsive, we observe a remarkably long lifetime of t = 700 ms 
(determined by an exponential fit). This lifetime is mainly limited by inelastic scattering 
of lattice photons. 
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Fig. 1. - Energy spectrum for a single particle in a ID lattice for three different potential depths 
Vo- The band width of the lowest Bloch band is given by 4 J where J is the hopping energy, a 
is the lattice period and E r = ir 2 h 2 /2ma 2 denotes the recoil energy. 




hold time (ms) 

Fig. 2. - left A state with two atoms located on the same site of an optical lattice has an 
energy offset « U with respect to states where the atoms are separated. Breaking up of the 
pair is suppressed due to the lattice band structure and energy conservation, so that the pair 
remains bound as a composite object, which can tunnel through the lattice. In the figure, n = 
denotes the lowest Bloch band and n = 1 the first excited band, right Breaking up of atom 
pairs in a shallow 3D optical lattice. [The potential depth is (10 ± 0.5) E r in one direction and 
(35 ± 1.5) E r in the perpendicular directions.] Shown is the remaining fraction of pairs for a 
scattering length of 100 ao (open diamonds) and a scattering length of about (0 ± 10) ao (filled 
circles) as a function of the hold time. The atom pairs quickly break up within a few ms if they 
do not interact, but stay together for a long time in case of repulsive interaction. The lines are 
fit curves of an exponential (dashed line) and the sum of two exponentials (solid line). 
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3. Analytical solution of two particle problem in an optical lattice 

3'1. General Discussion. - Our understanding of these stable pairs can be made more 
formal by an exact solution of the two particle Lippmann-Schwinger scattering equation 
on the lattice corresponding to the Bose-Hubbard Hamiltonian Eq. QJ. 

Denoting the primitive lattice vectors in each of the d dimensions by e^, we can 
write the position of the two atoms by x = 2i=i x i e i an d y = Sj=i 2/i e *7 where Xi,yi 
are integers, and we can write the two atom wave function in the form ^(x, y). The 
related Schrodinger equation from the Bose-Hubbard model [Eq. with homogeneous 
background, = then takes the form 



(2) 

where the operator 



-J ( A 



A, 



US 



x.y 



*(x,y) =£$(x,y), 



(3) 



A x *(x,y)=^ [*(x+e i ,y)+*(x-e i ,y) - 2tf(x,y)] 



denotes a discrete lattice Laplacian on a cubic lattice. Note that in order to express this 
in terms of the discrete lattice Laplacian we have added 4dJ4'(x, y) to each side of the 
Schrodinger equation. This effectively changes the zero of energy, so that E — > E + 4Jd. 
We then introduce relative coordinates r — x — y existing on the same lattice structure 
as the co-ordinate x, and center of mass coordinates R = (x + y)/2, existing on a lattice 
with the same symmetry as the original lattice but smaller lattice constant a/2. We then 
separate the wavefunction using the ansatz 



(4) 



*(x,y) = exp(iKR> K (r), 



with K the centre of mass quasi-momentum. This allows us to reduce the Schrodinger 
equation to a single particle problem in the relative coordinate, 



(5) 



2 J A* + Ek + US Tfi 



where Ej£ = 4J^ i= i [1 — cos(Ke»/2)] is the kinetic energy of the center of mass motion, 
and where the discrete lattice Laplacian for a square lattice is now given by 



(6) 



A^*(r) 



d 
i=l 



cos (Ke;/2) [*(r+ ei ) + *(r- ei ) - 2* (r)] 



The solutions of this Schrodinger equation can be found using the Greens function of 
the non-interacting problem with {7 = 0, which is defined by 



(7) 



[E-H }G K {E,r)=5 r ,o, 
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with S r o a three-dimensional Kronecker delta, and Ho = — 2 JA^ the Hamiltonian of the 
non-interacting system. This equation can be easily solved via Fourier transformation, 
G(E,r) = [l/(27r) d ] J d d kG(E, k) exp(ikr), and we obtain the solution 

(8) G K (E,k) 



E — eic(k) + ir\ 

where £k (k) accounts for the dispersion relation of the non-interacting system, 

(9) e K (k) = 4J^cos-^ [1 - cos^a)] . 

i=l 

The solutions of Eq. (J5J can be divided into two classes: scattering states, and bound 
(localised) states. We we will first analyze the scattering states. 

3'2. Scattering States. - Similarly to scattering problems involving particles in free 
space, the scattering states of particles on the lattice with energy E obey the Lippmann- 
Schwinger equation 

(10) fe(r) = i>%{v) + G K (E, r - r')V(r')^ E (r') 

r' 

with ip E — exp(zkr) an eigenstate of Hq with energy E = eic(k). In the present situation 
with a short range potential V(r) = US r ^o, the Lippmann-Schwinger equation can be 
solved via a resummation of the Born expansion and we obtain 

(11) Mr) = exp(ikr) - 8nJf E (K)G K (E, r) 
with scattering amplitude 

(12) fEiK) ~-^l-G K (E,0)U> 

where the total energy is E = ek,K + Ek, £k,K = 4J^i=i cos(Kej/2) [1 — cos(ke.;)], and 

4tt f dk 1 



(is) G K (£;,o) = gy 



2 J J {2n) d E /(2J) - 2 J2i=i cos (1 - cos ha) 

The scattering states ipE (r) correspond to two free atoms moving on the lattice and 
undergoing scattering processes. The corresponding energies appear as a continuum in 
Fig. |21 In order to make a connection to the scattering length in free space, we can 
consider the limit of small momenta of the incoming plane wave, i.e., k — > 0, K — > and 
E — > 0. Then the solution l|ll|l reduces in the limit r — > oo to 

(14) MrW o (r) + /(kjk0 exp^ 
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Fig. 3. - Two particle energy spectrum in a ID lattice for four different potential depths Vo as 
a function of centre of mass quasi-momentum K. The Bloch band for repulsively bound pairs 
is located above the continuum of unbound scattering states. The grey level for the shading of 
the continuum is proportional to the density of states. 



with the scattering amplitude 

1 U/(2J) 



(15) /(k,k') 



4tt 1 - otU/(2jy 



equivalent to the s-wave scattering length a s , while the constant a = lini£;_»o G(E, 0) is 
a re -0.25 Q3|. 

3'3. Bound States. - Note that the scattering amplitude in Eq. contains a pole, 
associated with a bound state. We now focus on these bound states in the regime, U > 0, 
which will correspond to a repulsively bound pair. First we note that we can write Eq. Q 
in the form 

(16) [E-H q ]G k (E,t) = — L 6 T ,oG K (E,r). 

As a consequence, the function ^^ s (r) = Gic(E,r) is a solution of the Schrodinger 
equation if the self-consistency relation is satisfied 

(17) U= ' 



Gk(E,0)' 
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which determines the energy Ebs of the bound state ip^. The resulting bound state 
wavefunction, ip^(r) ^ aus off exponentially for large r, and describes a bound two particle 
state travelling with center of mass momenta K through the lattice. The momentum 
distribution of this bound state is then given by Gk (Ebs j k) 

(18) VI s (k) = 



Ebs ~ £K(k) ' 

Note, that this wave function is not normalized. 

In three dimensions, Eq. i|17fl only has a solution for interaction strengths above a 
critical value, U > U c = — 2 J/G(0, 0) « 8 J, and thus we require U > U c for the formation 
of the bound two-particle state. The wavefunction ipjt(r) is then square- integrable, as 
shown in FigQJ For a deep lattice, i.e. U/J ^S> 1, bound pairs essentially consist of 
two atoms occupying the same site, whereas for small U/J, the pair is delocalized over 
several lattice sites. A main feature of the repulsive pair wavefunction is its oscillating 
character: the wavefunction amplitude alternates sign from one site to the next, as shown 
in Fig. 0J In quasimomentum space this corresponds to a wavefunction which is peaked 
at the edges of the first Brillouin zone, which is discussed in detail in section 

When motion is confined to one dimension the bound two particle state exists for 
arbitrarily small repulsive interaction U > 0, in contrast to the three-dimensional sit- 
uation. Here the energy of the bound pairs, computed from Eq. 1|17|) is Els(K) = 
2J[^/4(cos 4f) 2 + (U/2J) 2 + 2], which can be seen plotted in Fig. El as Bloch band of 
a stable composite object above the continuum of two particle scattering states. The 
figure shows how the binding energy (separation of these states from the continuum) 
increases as U/J is increased, and how the curvature of the band becomes less pro- 
nounced. In the limit of strong interaction, U 3> J the bound state energy reduces to 



8 



J. Hecker Denschlag A. J. Daley 



E(K) ~ 4 J + U + (4J 2 /i7)(l + cos Ka), which is consistent with our expectation of a 
positive binding energy U, and the center of mass energy of a composite object with an 
effective tunnelling matrix element J 2 /U . 

4. Numerical approach for repulsively bound pairs 

We want to be able to treat not just a single repulsively bound pair, but a lattice gas 
of many interacting repulsively bound pairs. This is important both in order to properly 
describe the effects of interactions on experimental measurements, and to investigate 
many-body effects on the behaviour of the pairs. Whilst perturbation theory can used 
to produce useful analytical models in some limits Jl] , it is possible in one dimension to 
treat the system in more general regimes using recently developed numerical techniques. 

4'f. Time- Dependent DMRG. - In one spatial dimension, the system of interact- 
ing repulsively bound pairs can be treated by directly simulating the Bose-Hubbard 
model time-dependently, using time-dependent Density Matrix Renormalisation Group 
(DMRG) methods |15l 1161 1171 ITS) . These methods allow near-exact integration of the 
Schrodinger equation for the ID many body system by an adaptive decimation of the 
Hilbert space, provided that the resulting states can be written efficiently in the form 
of a truncated matrix product state (this will be explained in more detail below) . This 
method has been successfully applied to several lattice and spin models of interest in con- 
densed matter physics, including systems realisable using cold atoms in optical lattices 
[ini EDI Ell 122 123 El ESI- The algorithm, both in the form originally proposed by Vidal 
and similar methods proposed by Verstrate and Cirac have also been generalised to the 
treatment of master equations for dissipative systems and systems at finite temperature 

HHE3 

4" 1.1. Matrix Product States. DMRG methods are based on a decomposition of 
the many-body wavefunction into a matrix product representation |18| . This requires 
that the original state can be expressed on a Hilbert space that is the product of local 
Hilbert spaces. For the Bose-Hubbard model, each local Hilbert space corresponds to 
a particular lattice site, and the basis states in the local Hilbert space correspond to 
different occupation numbers, running from to S — 1. We then write the coefficients of 
the wavefunction expanded in terms of local Hilbert spaces of dimension S, 

s 

(19) |*) = c ili2 ... tM \i 1 }(g>\i 2 )(S>...(E>\iM) 1 

i\ii...i M =^ 

as a matrix product state given by 

x 

(on) r . . . - "V rW ii\[i]r[ 2 l ia xPlrP] 12 rW]iM 

\* U J '-'l^-iM — / , i ai "ai 1 aia 2 rt a 2 1 a 3 a 4 ■■ - L a M -l ' 

Ol-»M-l 
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The r and A tensors are chosen so that the tensor \\J specifies the coefficients of the 
Schmidt decomposition P2JJ for the bipartite splitting of the system at site I, 

(2i) I^-EaEVL 1 -' 1 )!^ 1 -^ 1 ), 

a=l 

where xi is the Schmidt rank, and the sum over remaining tensors specify the Schmidt 
eigenstates, |</>q"'^) and |0q +1 "' M '). The key to the truncated matrix product state 
representation is that for many states corresponding to a low-energy in ID systems we 
find that the Schmidt coefficients A« , when ordered in decreasing magnitude, decay 
rapidly as a function of their index a |15) . Thus we can truncate the representation at 
relatively small x an d still provide an inner product of almost unity with the exact state 
of the system In implementations of this method we perform convergence tests for 
the state representation, that is, we vary the values of x an d S to check that the point 
at which the representation is being truncated does not affect the final results. We also 
make use of an optimised version of the code in which the Schmidt eigenstates are forced 
to correspond to fixed numbers of particles . This allows us to make use of the 

total number conservation in the Hamiltonian to substantially increase the speed of the 
code, and also improve the scaling with x an d S. With this number conserving code we 
are able to compute results with much higher values of x m a much shorter time than 
the original algorithm. 

4T.2. Time Dependence. Time dependence of these states can be computed for 
Hamiltonians acting only on neighbouring lattice sites because when an operator acts 
on the local Hilbert state of two neighbouring sites, the representation can be efficiently 
updated to the Matrix product state that best approximates the new state of the system. 
To do this, the T tensors corresponding to those two sites must be updated, a number 
of operations that scales as x 3 ^ 3 for sufficiently large x |15| . In this way, we represent 
the state on a systematically truncated Hilbert space, which changes adaptively as we 
perform operations on the state. The time evolution operator exp(— iHt), is then split 
into a product of operators, each of which acts only on a pair of neighbouring sites by 
means of a Suzuki- Trotter decomposition This is done in small timesteps St. 

Initial states can also be found using an imaginary time evolution, i.e., the repeated 
application of the operator exp(— HSt), together with renormalisation of the state. 

4"2. Numerical Investigation of Repulsively Bound Pairs. - In investigation of repul- 
sively bound pairs using these methods we are able to use parameters U and J directly 
corresponding to values of the lattice depth Vq in the experiments, and are also able to 
account for the background trapping potential. We typically study 10 — 30 pairs in 60 
lattice sites, and begin with an initial product state, corresponding to a random distribu- 
tion of doubly-occupied and unoccupied lattice sites. We then reduce the values of U and 
increase the value of J, using the same time dependence for the depth of the lattice Vb (t) 
as in the experiment. The single particle momentum distributions can then be calculated 
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efficiently from the matrix product state representation, and we also average the results 
over different initial configurations, to match the averaging over different ID tubes in the 
experiment. We can also perform lattice modulation spectroscopy, computing the time 
evolution of the many-body state when the parameters U and J vary as a function of 
time, based on the time dependence of the lattice depth V(t) used in the experiments 
(these calculations are similar to those in refs. |221 I25| ). 

5. Experimental realization 

In this section we describe the experimental steps to produce a pure ensemble of 
repulsively bound atom pairs in an optical lattice. We use 87 Rb as the atomic species 
for our experiments. 

5*1. BEC production. - We begin by creating a 87 Rb Bose-Einstein condensate (BEC) 
of 6 x 10 5 87 Rb atoms in spin state \F = 1, mp = —1) in a vacuum apparatus featuring a 
magnetic transport line |30U31| . This transport initially transfers laser cooled atoms from 
the chamber of the magneto-optical trap (MOT) into a UHV glass cell (pressure < 10~ n 
mbar) which offers good optical access from all sides. Here the BEC is produced in a 
QUIC trap [3211211 wu - n trapping frequencies 0J x , y ,z = 2n x (15, 15, 150) Hz at a magnetic 
bias field of 2 G. Afterwards the QUIC trap is converted into a Ioffe-type magnetic trap 
with trap frequencies to XtVtZ — 2tt x (7, 19, 20) Hz) by adjusting the currents through 
the quadrupole and Ioffe coils and by applying additional magnetic field gradients. This 
moves the BEC over a distance of 8 mm into the center of the QUIC quadrupole coils 
which are later used to generate the homogeneous magnetic field for Feshbach ramping. 

5'2. Loading into lattice. - Within 100ms the BEC is adiabatically loaded into the 
vibrational ground state of an optical lattice which is 35 E r deep (E r = 2n 2 h 2 /m\ 2 , 
where m is the mass of the atoms) . Our 3D lattice is cubic and consists of three pairs 
of retro-reflected intensity-stabilized laser beams, which propagate orthogonally to each 
other. They are derived from a frequency-stable single-mode Ti: Sapphire laser (< 500 
kHz linewidth) with a wavelength of A = 830.44 nm. For this wavelength, the laser 
is detuned by about 100 GHz from the closest transition to an excited molecular level, 
minimizing light induced losses as a precondition for long lifetimes of pairs and molecules. 
The laser beams are polarized perpendicularly to each other and their frequencies differ 
by several tens of MHz to avoid disturbing interference effects. The waists of all three 
beams are about 160 /im and the maximum obtainable power is about HOmW per beam. 
At this stage about 20% of the condensate atoms are grouped in pairs of two into the 
lattice sites. 60% of the condensate atoms are found in singly occupied sites, and another 
20% percent of atoms are located in triply and more highly occupied lattice sites 33 . 

5'3. Purification scheme. - In order to remove all atoms from those lattice sites that 
are not occupied by exactly two atoms we use a purification scheme which involves an 
intermediate step in which Feshbach molecules are produced. For this, we turn off the 
magnetic trap and flip the spins of the 87 Rb atoms from their initial state \F — 1, mp — 
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— 1) to \F — 1, mp — +1) by suddenly reversing the bias magnetic field of a few G. This 
spin state features a 210 mG wide Feshbach resonance at 1007.40 G |33|. By adiabatically 
ramping over this resonance we convert pairs of atoms in multiply occupied lattice sites 
into Rb2 Feshbach molecules with almost unit efficiency [33]. After the Feshbach ramp, 
fast inelastic collisions will occur in sites that were initially occupied with more than 
two atoms. These exothermic collisions between either a created molecule and an atom 
or between two molecules will remove these particles from the lattice. At this stage the 
lattice consists only of sites which are either empty, filled with a single atom, or filled with 
a single Feshbach molecule. A subsequent 3ms long combined radio-frequency (rf) and 
optical purification pulse removes all chemically unbound atoms, thus creating a pure 
molecular sample of about 2 x 10 4 molecules (231 - The microwave drives the transition 
at a frequency of 9113MHz between levels which correlate with \F — I, nip — +1) and 
\F = 2. nip = +2). The light pulse drives the closed transition \F = 2, mp = +2) — > 
\F = 3, nip = +3). The optical transition frequency is 1402 MHz blue detuned compared 
to the transition at zero magnetic field. The light literally blows the atoms out of the 
lattice, while the direct effect of the microwave and light field pulse on the molecules 
is negligible because the radiation is off resonance. Finally, sweeping back across the 
Feshbach resonance we adiabatically dissociate the dimers and obtain a lattice where 
2 x 10 4 sites are filled with exactly two atoms. According to section 131 at the deep lattice 
depth of 35 E r (U/J « 3700), the corresponding two atom wavepacket matches perfectly 
with the wavefunction of the repulsively bound atom pair. By adiabatically lowering the 
lattice depth (typically within a few ms) in a horizontal direction we can then produce 
ID repulsively bound atom pairs states for arbitrary values of U /J. 

6. Experiments 

In the following we discuss the properties of the repulsively bound pairs that were 
experimentally investigated by measuring their lifetime, quasi-momcntum distribution 
and binding energies. By varying the effective interaction between the atoms with the 
help of the Feshbach resonance we can also create lattice induced bound atom pairs which 
are based on attractive interactions. 

6T. Pair lifetime. - We have already seen from Fig. fright) that for a repulsive 
interaction with a s — 100 ao the lifetime of the pairs is remarkably long (700 ms, expo- 
nential fit). This lifetime is mainly limited by inelastic scattering of lattice photons |33|. 
and greatly exceeds the calculated time for an atom to tunnel from one site to the next, 
2itTi/{AJ) ~ 4ms. The lifetime measurements are based on lowering the lattice depth to 
a chosen height, and then measuring the number of remaining pairs after a variable hold 
time. In order to make this measurement, the lattice is adiabatically raised again to its 
full initial depth of Vq = 35 E r . Using the Feshbach resonance atoms in doubly occu- 
pied sites are converted to Feshbach molecules with near unit efficiency |33) . and another 
combined rf-light purification pulse then removes all remaining atoms (which stem from 
now dissociated pairs) as in the original preparation step. Afterwards the molecules are 
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again converted back into atoms, and can then be detected via conventional absorption 
imaging. 

6'2. Quasi-momentum distribution. - We have experimentally investigated the quasi- 
momentum distribution of the pairs in various regimes by mapping it onto a spatial 
distribution, which we measured using standard absorption imaging. For this, we first 
adiabatically lower the lattice depth in a horizontal direction at a rate of 1.3 E r /ms 
to a pre-chosen height while the lattice depth in the other two directions are kept high 
(35 E r ). This will prepare repulsively bound pairs at the chosen lattice depth. We 
then turn off the lattice rapidly enough so that the quasi-momentum distribution cannot 
change, but slowly with respect to the bandgap, so that single-particle quasi-momenta 
are mapped to real momenta |31II35| . We have typically employed linear ramps with rates 
of 0.2 E r /fis. The resulting momentum distribution is converted to a spatial distribution 
after ~ 15 ms time of flight. Fig. [S] shows two measured quasi-momentum distributions 
for lattice depths Vq = 6(20), respectively. The top row shows the bare images of the 
atomic density taken in the laboratory. Below are the corresponding quasi-momentum 
distributions in horizontal direction. It is clearly visible that for low lattice depths the 
quasi-momentum distributions are peaked at the edges of the first Brillouin zone. For 
deep lattices, however, the first Brillouin zone is homogenously filled and the quasi- 
momentum distribution has a flat top shape. This latter distribution is reminiscent of 
the one observed for a dephased ensemble of ultracold atoms in the lowest band of a 
lattice The agreement between experimental data and theoretical calculations is 

quite good. Note, however, that the experimental distributions appear to extend beyond 
the first Brillouin zone. This is an experimental artifact related to repulsion between 
atoms during expansion (before imaging) and also relatively long imaging times (many 
photons are scattered from each atom, which performs a random walk). This leads to 
smearing out of the sharp structure at the edge of the Brillouin zone. Fig. shows in 
a more continuous fashion the dependence on lattice depth Vo of the quasi-momentum 
distribution for repulsively bound pairs for both experiment and numerical simulation. 
As discussed before, the peak structure is more pronounced for lower values of Vo, and 
diminishes for larger Vo- This characteristic is a clear signature of the pair wavefunction 
for repulsively bound pairs. 

It is important to note that in all cases here we measure the distribution of single- 
atom quasi-momenta in a large sample. That we still obtain the peaked distribution 
characteristic of repulsively bound pairs is non-trivial. In fact, if we just take a single 
repulsively bound pair with centre of mass quasi-momentum K ^ 0, its single- atom 
momentum distribution will not be peaked anymore at the edges of the first Brillouin 
zone. The peak will be somewhat translated towards the center of the first Brillouin 
zone. Fortunately, with increasing \K\, the peak in the single-particle quasi-momentum 
distribution also becomes less pronounced. As a result, when we average over a roughly 
uniform distribution of centre of mass quasimomenta K for a dilute gas of repulsively 
bound pairs, we still observe the pronounced peaks at the edges of the Brillouin zone. This 
is confirmed by the numerical simulations and their good agreement with experiments. 
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Fig. 5. - Quasi-momentum distributions of atoms in the lattice for (a) Vo = 5Er and (b) 
Vo = 20Er. above Images show absorption images of the atomic distribution after release 
from the 3D lattice and a subsequent 15 ms time of flight. The horizontal and vertical dashed 
lines enclose the first Brillouin zone, below Corresponding quasi-momentum distributions in 
the horizontal-direction, after integration over the vertical-direction. For comparison numerical 
simulations (see text) are also shown (dashed lines). The density values have been scaled to 
facilitate comparison between experimental and theoretical results. 
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Fig. 6. - Momentum distributions similar as the ones shown in Fig|S| are plotted here as a 
function of lattice depth Vo- Left: experiment. Right: numerical calculation. 
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Fig. 7. - Modulation spectroscopy of repulsively bound pairs, (left) By modulating the optical 
lattice amplitude with the proper frequency, the pair can dump its binding energy into the 
lattice motion and subsequently break up. (right) Typical resonance dip showing the remaining 
number of atom pairs as a function of the modulation frequency, for Vo ~6E r . The black line 
is a Gaussian fit, a choice which was justified by numerical calculations. 



2.5 




Fig. 8. - Measured resonance frequencies of the modulation spectroscopy as a function of the 
lattice depth. The resonance frequency was determined from resonance curves similar to the one 
in Fig. (right). The experimental points (filled circles) show good agreement with numerical 
simulations (crosses) and also coincide with the onsite collisional energy shift U (line). Experi- 
mental error bars correspond to the 95% confidence interval for the Gaussian fit parameters of 
the resonance dips. 
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6'3. Modulation spectroscopy . - By modulating the depth of the lattice at a chosen 
frequency we can determine the binding energy of the pairs (see Fig. left). For 
appropriate modulation frequencies, the pairs can dump their binding energy into the 
lattice motion and dissociate. Fig.0(right) shows a typical resonance curve of the number 
of remaining pairs as a function of the modulation frequency. The resonance frequency 
of about 1.5 kHz (for a lattice depth of Vq — 6 E r ) agrees well with the calculated 
binding energy of a pair. The width of the resonance curve can be understood, as the 
pair will decay into a continuum of scattering states which has an energy width of up to 
8 J [depending on the initial centre of mass quasimomentum K (see Fig. 15)]. In addition 
to this width, broadening due to Fourier limited modulation pulses and inhomogeneity 
effects in the lattice will occur. 

Modulation spectroscopy measurements were carried out for a variety of lattice depths 
(see Fig. |SJ). The resonance positions are in good agreement with numerical simulations 
and essentially coincide with interaction energy, U. 




B(G) a(ao) 

1007.5 (-110) 
1007.6 (0) 
1007.7 (+30) 
1007.8 (+50) 



0.2 0.4 0.6 
B- 1007 G (G) 



Quasi-momentum (jt/a) 



Fig. 9. - From repulsively to attractively bound atom pairs, left With the help of a Fesh- 
bach resonance around 1007.4 G, we can choose the effective interaction of the paired atoms 
by controlling the scattering length a. ao is the Bohr radius. The inserted images show quasi- 
momentum distributions similar to the ones of Fig. For effectively attractive interaction the 
quasi-momentum distribution is peaked around momentum, right The momentum distribu- 
tion for atom pairs as a function of magnetic field (scattering length.). At zero scattering length 
the distribution has a flat top shape. The shown data correspond to experiments where the 
lattice depth Vq had been adiabatically lowered in ID below 3E r . 



6'4. Attractively bound pairs. - Making use of the Feshbach resonance at 1007.40 G 
we can tune the effective interaction of the atoms within the pair (see Fig. HO left). It is 
then possible to also create bound atom pairs which are based on attractive interaction. 
After initial production of repulsively bound atom pairs in the deep lattice (Vq — 35E r ), 
we applied an appropriate nearly homogenous magnetic offset field. This tuned the 
scattering length of the atomic pair from its default value of a s = lOOag to negative 
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scattering length of up to a s = — llOao- Afterwards the optical lattice was lowered 
as before. In contrast to repulsively bound pairs where the momentum distribution is 
peaked at the edges of the first Brillouin zone, the momentum distribution for attractively 
bound pairs is peaked in the center of the first Brillouin zone. This goes along with 
having a bound state with minimal internal energy. Fig. (right) shows how the quasi- 
momentum distribution of the pairs changes continuously as the scattering length is 
changed. Interestingly, for non-interacting atoms (a s = 0) the distribution again becomes 
a flat top shape. 

With respect to stability, we find that lifetimes of bound atom pairs are similar for 
scattering lengths of equal size but opposite sign. 

7. Repulsively bound pairs of fermions 

Although so far we have only discussed repulsively bound pairs which are composed 
of bosonic atoms, analogous pairs can also be formed from fermions or boson-fermion 
mixtures. These systems will exhibit interesting physics based on their composite nature 
and the quantum statistics of their components. 




Fig. 10. - Possible realization of repulsively bound pairs with fermionic atoms, left In the case 
of two identical fermions (same spin) the atoms have to be in different bands due to the Pauli 
exclusion principle, right Fermions of different spins or Bose-Fermi mixtures, however, can 
occupy the same band. 



Of course, in a single species fermion experiment it is not possible to put two identical 
fermions into the same site and band due to the Pauli exclusion principle. The two 
atoms would have to be at least in different bands (see Fig. left) , and even then the 
interaction between them typically would be very small in the ultracold regime. 

These problems do not arise using a two-component spin mix of fermions (see Fig. 
1101 right), as two atoms of different spin can share the same site and band and can also 
interact strongly. The fact that a higher site occupancy than two is again strictly forbid- 
den could be advantageous in the initial production of pairs. Furthermore, the pairing of 
two fermions can result in a pair with bosonic character. In a 3D environment pairing of 
fermions recently lead to interesting experiments studying the BEC-BCS transition (see 
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e.g. [3S1I2ZH3H1 EH! EHB HZ])- It would be interesting to study similar properties to this 
transition with repulsive pairing, investigating the system as the interaction strength is 
changed. 

Also it would be very interesting to study repulsive pairs which are a composite objects 
of a fermion and a boson. One question would be how these bound states, would interact 
with each other, and how, for example, the bosonic atoms within the pair would mediate 
next neighbor interactions 

8. Other related physical systems 

Although no stable repulsively bound pairs have previously been observed, their 
physics is partially related to other physical systems. Here we briefly discuss a few ex- 
amples of such systems with bound states which are in a way reminiscent of repulsively 
bound states. 

8'1. Pairing resonances in many-body systems. - For example, resonance behaviour 
based on similar pairing of Fermions of different spin in the Hubbard model was first 
discussed by Yang 021 , and plays an important role in SO(5) theories of superconductivity 
|4*5| . There are several examples of many-body bound states that can occur for repulsive 
as well as attractive interactions, such as the resonances discussed in the context of the 
Hubbard model by Demler et al. |44|. Such resonance behaviour is common in many- 
body physics, although states of this type are normally very short-lived. Optical lattice 
experiments will now provide an opportunity to prepare and investigate stable versions 
of such states, which until now have only appeared virtually as part of complex processes. 

8'2. Excitons. - The stability and many-body physics of repulsively bound pairs is 
perhaps most closely associated with that of excitons, which are attractively bound pairs 
of a particle in the conduction band and a hole in the valence band of a periodic system 
|45|. These bind to form a composite boson, a gas of which can, in principle, Bose- 
condense. Excitons are excited states of the many-body system, but are bound by an 
attractive interaction between the particle and hole that form the pair. They are also 
discussed in the specific context of fermionic systems. However, a single exciton on a 
lattice could have a description very similar to that of a single repulsively bound pair, 
and could be realised and probed in optical lattices experiments 46 . 

8"3. Photonic crystals. - Repulsively bound atom pairs in an optical lattice are also 
reminiscent of photons being trapped by impurities in photonic crystals |47| , which con- 
sist of transparent material with periodically changing index of refraction. An impurity 
in that crystal in the form of a local region with a different index of refraction can then 
give rise to a localized field eigenmode. In an analogous sense, each atom in a repulsively 
bound pair could be seen as an impurity that "traps" the other atom. 

8'4. Gap solitons. - An analogy can also be drawn between repulsively bound atom 
pairs and gap solitons, especially as found in atomic gases 1491 1501 15 1| . Solitons 
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are normally a non-linear wave phenomenon, and in this sense have a very different 
behaviour to repulsively bound pairs, which exhibit properties characteristic of many- 
body quantum systems. However, there has been increasing recent interest in discussing 
the limit of solitons in atomic systems where very few atoms are present, giving rise to 
objects that are often referred to as quantum solitons |52l 1581 151] . These are N-body 
bound states in ID, and thus a 2-atom bright quantum soliton is a bound state of two 
atoms moving in ID. In this sense, the solution for a single repulsively bound pair in ID 
is related to a single quantum soliton on a lattice. 

9. Conclusion 

We have reviewed theoretically and experimentally the physics of repulsively bound 
pairs of atoms in an optical lattice. The good agreement between experiment and the- 
ory exemplifies the strong correspondence between the optical lattice physics of ultra- 
cold atoms and the Hubbard model on a new level, a connection which has particular 
importance for applications of these cold atom systems to more general simulation of 
condensed matter models and to quantum computing. The existence of such metastablc 
bound objects will be ubiquitous in cold atoms lattice physics, giving rise to new potential 
composite objects also in fermions or in systems with mixed Bose-Fermi statistics. These 
states could also be formed with more than two particles, or as bound states of existing 
composite particles. Repulsively bound pairs have no direct counterpart in condensed 
matter physics due to the strong inelastic decay channels observed in solid state lattices, 
and could be a building block of yet unstudied quantum many body states or phases. 
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